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Motivated by recent experiments in Bose-Einstein condensed atoms that have been confined in 
toroidal traps, we examine the stability of persistent currents in such systems. We investigate the 
extent that the stability of these currents may be tunable, and the possible difficulties in their 
creation and detection. 
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I. INTRODUCTION 

The recent advances in the physics of cold atoms have 
allowed experimentalists to engineer many of the prop- 
erties of these systems. Remarkably, the experiments of 
Refs. [H, 01 have managed to trap atoms in toroidal traps, 
while in Ref. Q persistent currents were created and ob- 
served in toroidal traps. Furthermore, in an older exper- 
iment, Ref. [3] has investigated the dissipationless flow 
of an obstacle in an elongated Bose-Einstein condensate 
below some critical velocity. 

Such simple trapping geometries makes these systems 
very appealing mainly for two reasons, one theoretical 
and one technological. Theoretically these gases are ideal 
for testing fundamental superfluid properties (such as 
persistent currents, for example), and they may real- 
ize exactly-solvable one-dimensional models. Technologi- 
cally, the possibility of creating persistent currents whose 
stability is tunable externally may lead to important, or 
even revolutionary applications. 

Motivated by these facts, we focus in the present study 
on the stability of persistent currents in a quasi-one di- 
mensional, toroidal trap. Three are the classes of prob- 
lems that we consider. In the first one, we consider the 
situation where one manipulates the trap appropriately, 
in order to achieve the desired effect on the stability of 
currents. In the second, we investigate unavoidable com- 
plications in realistic experiments, which may affect the 
stability of persistent currents. The third class includes 
the problems associated with the detection of persistent 
currents. 

Many aspects of the problem that we examine in 
our study have been investigated thoroughly in previ- 
ous studies. We do not attempt to give a complete list 
of references, but rather we mention just few of them 
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In what follows, we first consider our model in Sec. II. 
Then, in Sec. Ill we consider a step potential along the 
torus and investigate the conditions that destabilize the 
current most easily. In Sec. IV we consider a periodic po- 
tential along the torus and examine the stability of per- 
sistent currents in such a potential. Then, in Sec. V we 



consider a random potential and via a statistical analysis 
of our results, we examine the effect of the length scale 
of the irregularities on the stability of the current. In 
Sec. VI, we investigate the effect of gravity on a tilted 
torus, and in Sec. VII we propose a method which al- 
lows the detection of a current. Finally, in Sec. VIII we 
present the conclusions of our study. 



II. MODEL 

In our model we consider one-dimensional motion of 
the atoms with periodic boundary conditions, which cor- 
responds physically to the motion of the atoms in a very 
tight toroidal trap. Quasi-one dimensional motion is ac- 
complished as long as the interaction energy is much 
smaller than the excitation energy transversely to the 
axis of the torus. We also assume for the atom-atom colli- 
sions the usual contact potential, Vint(r— r') — UoS{r~r') 
with Uo = inh^ / M . Here, age is the scattering length 
for elastic atom-atom collisions and M is the atomic 
mass. 

The mean-field approximation that we use in this study 
implies that the condensate order parameter '^(0) satis- 
fies the nonlinear, Gross-Pitaevskii-like equation [20l[2]|. 

-^ + v{e)^ie) + 27rj\^f{e)\H^tiM>, (i) 

where we have set fi = 2M — R — 1, with R being 
the radius of the torus. Here 6 is the azimuthal an- 
gle, V{9) is the external potential, and /i is the chem- 
ical potential. Both V{9) and fi are measured in units of 
the kinetic energy T = /{2MR^). The ratio between 
the interaction energy and the kinetic energy is equal to 
7 = naUa/T = ma^^R/S. Here tiq = N/{2nRS) is the 
average atom density, iV 3> 1 is the atom number, and 
S is the cross section of the torus (with R ^ VS). As 
long as 7 ^ iV^, the system is away from the Tonks- 
Girardeau limit of impenetrable bosons [13], that we do 
not consider in our study. 

To investigate the stability of persistent currents, we 
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use the time-dependent version of Eq. ([T]) 

^^^--^ + Vm{0)+2nj\^{e)\'^. (2) 

We propagate some initial state in imaginary time, mak- 
ing the substitution t = it 23], and solve the equation 

- = -0 + Vi9)^i9) + 2TT^\^i9)\H - fi^. (3) 

More specifically, in the above equation we specify 
^'(6', T — 0) and look for the convergent solution that 
emerges for long enough times from the above time- 
dependent equation. 

The states that correspond to successive values of (the 
quantized) circulation may be separated by an energy 
barrier, provided that the interaction is strong enough 
[23 , [U, [2^ . These barriers give rise to stable persistent 
currents. Therefore, if one starts with the initial condi- 
tion '5(6', T = 0) = e'^/\/27r (for example) that has one 
unit of circulation and a uniform density distribution, if 
there exists an energy barrier between this state and the 
current-free state, then the long-time solution of Eq. ([3]), 
still has one unit of circulation, with possibly some angle- 
dependent variation in its density, depending on the form 
of the potential V{9). On the other hand, in the absence 
of an energy barrier between the two states (which is 
the case for sufficiently weak, or attractive interactions) 
the system is energetically unstable, and it relaxes to the 
circulation-free state. 

The criterion for the stability of circulation that wc 
use in our study is equivalent to a condition that resem- 
bles the Landau criterion for superfluid flow More 
precisely, in order to have stability of the superflow, the 
"drift" velocity of the atoms u must not exceed the speed 
of sound c. For example, for one unit of circulation that 
we consider here, the drift velocity is equal to h/{MR). 
Also, in the limit of weak interactions the speed of sound 
is given by j26|] (see, e.g., Eq. or Ref. f27||), 

c=— ^(1 + 7), (4) 

and in order for u to be smaller than c, 7 has to be 
larger than unity, in agreement with the criterion of en- 
ergetic stability. Finally, we mention that while we have 
restricted our study to the transition between states of 
circulation equal to unity and zero, similar effects occur 
in transitions between states of higher, successive values 
of the circulation. 



III. DESTABILIZATION OF PERSISTENT 
CURRENTS WITH A STEP POTENTIAL 

As a first application of our method, we consider a 
step-like potential of a fixed height (equal to unity in our 
units) along the torus. 



and investigate the critical value 7c of the coupling 7 that 
is necessary to maintain the stability of the current, as 
function of 9o. The form of the potential that we choose 
is not accidental. Among all the possible functional forms 
of V{9) that we investigated with V^nax ~ Knin — 1, this 
is the one that requires the highest value of 7c, for a 
static potential. Therefore under these conditions, at 
least within the functional forms of the potentials that 
we have considered, the derived value of 7c may serve as 
an upper bound for the critical coupling that is neces- 
sary to sustain a persistent current of circulation equal 
to unity. 

Following the method that was described in the previ- 
ous section, we plot in Fig. 1 7c versus 9o. As one sees 
in this plot, the highest value of 7c that is necessary to 
stabilize the current occurs for 9q/{2tt) ~ —0.18. The 
position of 9^ decreases with increasing Knax — Knin- 

For ^0 ~* "I" J E^nd ~^ 7c tends to 3/2, which 
is the critical value of 7c that corresponds to a constant 
V[9) ^lUll^l. For intermediate values of 6*0 the situa- 
tion becomes more interesting, as in this case the length 
scale of variation of V{9) is comparable to the coherence 
length ^ and as a result 7c is higher than 3/2. More pre- 
cisely, when 7 is of order unity, the coherence length ^, is 
~ iJ, since ^/i? = 7^^/^. The fact that ^ is comparable 
with the radius of the torus, and thus comparable with 
the length scale of variation of V{9), implies that the cur- 
rent becomes fragile. A crude guess for 9q/{2tt) that re- 
quires the maximum value of 7c is the middle of the torus 
(i.e., when the step potential extends over one half of the 
torus), as we also discuss in the toy model that follows 
below. We mentioned earher that for T4iax — K„in = 1 the 
actual value is ~ —0.18, with the corresponding value of 
7c being ~ 3.8. 

A simple toy model for this problem gives a qualita- 
tively correct answer for 7c as function of ^o- Let us 
consider just the states $o(^) = l/\/27r and $i((?) = 
e*^/\/27r in the order parameter, in the form 

*triai(0) = \/T~?3'o(f?)+e^M$i(0), (6) 

where A is some phase factor. The value of A is de- 
termined from the minimization of the potential energy, 
which turns out to be A = —{9q + tt)/2. Clearly the above 
state has an expectation value of the angular momentum 
per particle equal to I. The corresponding density is 

ntriai(e) = 7^[1 + 2V^(l-Osin(6^ - 9ol2)l (7) 
zvr 

and also the (minimized) energy per particle is, for the 
step potential of Eq. ([5]) , 

£;-2-^ = (l+7)/-7^^-yi(r37)^cos(|). (8) 

From the above dispersion relation we find the critical 
value 7c that gives rise to a local minimum close to Z — 1 . 
The specific (highly truncated) order parameter gives 
7c = 1 for a constant potential (i.e., for = i""); while 



71 



1 2 3 4 5 6 
m 



FIG. 1: The critical value of the coupling constant 7c versus 
the location of the step So, in a step potential of the form of 
V{0) defined in Eq. ([5|. The value of 7c = 3/2 for the uniform 
case (dashed line) is also shown as a reference. 

7c has a maximum for 6*0 = 0, equal to w 2.295. In addi- 
tion, 7c(^o) is symmetric around — 0. In the limit of 
weak disorder and weak coupling, where the actual order 
parameter is closer to Atrial (6*), the density profile of the 
cloud is sinusoidal, and the step potential just determines 
the position of the maximum/minimum of the density, as 
one can see in Eq. ([7]). As a result, 7c(— 6*0) = Ici^o)- For 
higher couplings and higher values of Vmax ~ Knin, the 
density of the cloud gets distorted from the sinusoidal 
form, as it localizes more in the region of lower potential 
energy. The function 7c (^0) becomes then asymmetric 
with respect to 9o = 0, as shown in Fig. 1 (which shows 
the full numerical calculation.) 

IV. STABILITY OF PERSISTENT CURRENTS 
IN THE PRESENCE OF A PERIODIC 
POTENTIAL 

Another interesting question is the stability of persis- 
tent currents in the presence of a periodic potential that 
acts along the torus. Experiments with periodic poten- 
tials created by optical lattices have already been per- 
formed [2^ . Numerous theoretical studies have examined 
this problem, too. This is a very interesting problem, 
since the presence of a periodic external potential, com- 
bined with the effect of the interactions give rise to novel 
states. For example, we refer to Refs. [2^, HO, HH, [H, 
for studies of the Gross-Pitaevskii equation in the pres- 
ence of a periodic potential, for calculations of the band 
structure, and for the study of the superfluid properties 
of these systems. 

In our study we consider a sinusoidal external potential 
of the form 

V{9) = cos(m6i), (9) 

where m = 1,2,3,..., as the condition V{d + 2-k) — V{9) 
requires. This potential has a period equal to 27r/m. 
Again, we examine the stability of a current with one unit 



FIG. 2: The critical value of 7c versus the parameter m (dots), 
for the periodic potential of Eq. @. The solid line is a guide 
for the eye. The value of 7c = 3/2 for the uniform case 
(dashed line) is also shown as a reference. 

of circulation, starting with ^'(6', t = 0) = e*^/\/27r and 
finding the critical value of 7 that gives rise to a stable 
current, for various values of ni — 1, ... ,6. The results 
of this calculation are shown in Fig. 2. As m increases, 
the value of 7c that is required to give rise to a stable 
current decreases. 

For m = 1, the potential is qualitatively similar to a 
step potential examined in Sec. Ill, but it is smoother 
(which decreases 7c) and has a larger Vmax — — 2 
(which increases 7c). For a step potential with T^nax — 
T^nin = 2, that is located in the middle of the torus 
(6'o = 0), one finds 7c — 4.8, in rough agreement with 
the maximum value of 7c shown in Fig. 2. Furthermore, 
as m increases, eventually the potential becomes homo- 
geneous (on the length scale of variation of the den- 
sity), in which case 7c = 3/2 [13, IHl, 113]. This is 
clearly seen in Fig. 2. The crossover region to the ho- 
mogeneous case is accomplished when the "lattice con- 
stant" 7r/m, associated with the periodic potential V{6), 

IS ^ S,/ R = ■ The above equation implies that 

TO TTjl^^ = 7r(3/2)^/^ w 3.9, in (rough) agreement 
with our result shown in Fig. 2. 

V. STABILIZATION OF A CURRENT IN THE 
PRESENCE OF A RANDOM POTENTIAL 

In a previous study [l^l , we investigated the value of 7c 
that is necessary to stabilize a current in the presence of 
a piecewise constant potential, with a randomly chosen 
amplitude. Here we present data that we have collected 
from a statistical analysis of our simulations, where we 
examine 7c as function of the number of steps that we 
choose in V{9). 

More specifically, we compute the average value of 7c, 
as well as the standard deviation cr(7c), for 1000 different 
random step potentials, for which there are s steps, with 
a width (27r/s). The value of the potential within each 
step is drawn uniformly from the interval [— 1 , 1] . Figure 
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FIG. 3: The average value of the critical coupling (7c) (dots), 
as well as the standard deviation (bars), versus the number 
of steps s in a randomly-chosen potential. These results are 
derived from 1000 different random potentials for each num- 
ber of steps s = 5, 10, 15, 20. The value of 7c = 3/2 for the 
uniform case (dashed line) is also shown as a reference. 

3 shows the result of these calculations. We observe that 
(7c) is a decreasing function of s. 

As in the case of a step potential, one may argue that 
as s increases, the length scale of variation of the ran- 
dom potential becomes increasingly smaller than the co- 
herence length, and eventually one recovers the result 
7c = 3/2 of the homogeneous torus, when V{9) is con- 
stant. On the other hand, as s decreases, eventually the 
two length scales become comparable and the random 
potential destabilizes the current more easily, requiring a 
higher value of 7c for the current to become stable. 
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FIG. 4: The critical value of 7c versus log^g{5), for R — 1 mm 
(left), 0.1 mm (middle), and 0.01 mm (right). The value of 
7c = 3/2 for the uniform case (dashed line) is also shown as 
a reference. 

(which correspond to the three different values of R) , are 
almost identical, as they can be obtained approximately 
by a shift in 6, which is cx R^^. 

It is instructive to get an estimate for the energy scale 
MgR in Eq. (jlOp . Considering Rubidium atoms, for ex- 
ample, and for a radius of the torus R = 0.1 mm, MgR is 
on the order of /LtK, which is a rather large energy scale. 
Therefore, in order for the effect of gravity to be unim- 
portant, the torus has to be tilted slightly, in which case 
sin (5 < 1. 



VII. DETECTION OF PERSISTENT 
CURRENTS 



VI. EFFECT OF GRAVITY ON THE 
STABILITY OF PERSISTENT CURRENTS 

One question that needs to be investigated within the 
problems that we examine here is the effect of gravity 
on the stability of persistent currents. Already in other 
experiments with cold atoms, gravity has played an im- 
portant role, see e.g. [l^]. In real life, the torus cannot 
be perfectly horizontal, and as a result, there is an angle- 
dependent potential that acts on the atoms due to the 
gravitational force. This potential is 

V{e) = [MgR sin S) cos 9, (10) 

where g is the acceleration of gravity and S is the angle 
between the plane of the torus and the horizontal plane. 

Figure 4 shows 7c as function of \ogiQ{d), for three 
values of i? = 1 mm (left), 0.1 mm (middle), and 0.01 
mm (right). As the angle S increases for a given R, 7c 
increases too, since the effect of the gravitational field 
gets more important. For the same reason, for a fixed 
angle 5, as the radius of the torus R increases, the critical 
value of 7c increases, too. For the small values of the 
angle S that we have considered, sin 5 « i5, and therefore 
the shape of the three different curves that are shown 



A serious issue in these experiments is to be able to 
know whether there is circulation in the gas, or not. This 
question can be resolved with use of interference tech- 
niques (ssj . however it would be convenient to have also 
a way to measure the circulation using other methods. 

It is natural to think of the single-particle density 
distribution n{d) = |^'(0)p as a possible way to deter- 
mine the value of the circulation. Any constant potential 
V{9) results in a homogeneous density distribution, as 
the order parameter is e"'^/\/27r, for any value of the 
circulation n{2TTh/M). On the other hand, if one uses 
a probe potential V{6) that is spatially-dependent, the 
density distribution n{d) does depend on the value of k. 
Therefore, measuring the density in the presence of some 
spatially-dependent potential V{6) may allow us to de- 
termine the value of the circulation. 

We thus consider such a probe potential in the form 
of a Gaussian unity dimple (which may be realized via 
some laser beam) with (rms) width Wd 

F(0) = -exp(-(?V2«;2). (n) 

Figure 5 shows the density distribution n{6) in the pres- 
ence of V{9) of Eq. PT|) . with Wd = 0.2 for the states with 
zero circulation (k — 0), and unit circulation [k, = 1). 
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FIG. 5: The (normalized) density distribution n(Q) of the 
state of lowest energy, in the presence of ViQ) given by 
Eq. (|lip . for zero circulation, k = (dashed curve), and 
unit circulation k = \ (solid curve), for a coupling strength 
7 = 2.2. The horizontal dotted line shows the density of the 
uniform state, n(Q) — 1/2tt. All densities are measured in 
units oi N/{RS). 

The chosen value of the coupling is larger than the crit- 
ical one, 7 = 2.2 > 7c ~ 2.1. The state with k = 1 has 
a higher kinetic energy, and as a result the classically- 
forbidden region is more narrow, resulting in a wider 
density distribution around the minimum of the probe 
potential V{6). Ideally one would like the density pro- 
files corresponding to different values of the circulation 
to be as different as possible, even in the limit of a very 
weak probe potential |V^(6')| ^ 1. 

To quantify the difference between the two density dis- 
tributions illustrated in Fig. 5, one may introduce the 
difference between the highest and the lowest value of 
the density in each case. An = rimax — 'T-min = n{9 = 
0) — n{9 = ±7r). Remarkably, this quantity differs by 
more than a factor of two (« 2.7) between the two den- 
sity profiles that we consider in the specific example. 



The above results indicate that a possible way to mea- 
sure the circulation (non-destructively), would be to turn 
on the probe potential adiabatically, measure the density 
and extract the value of n (non-destructively) , and finally 
turn off the probe adiabatically again. 



VIII. SUMMARY AND CONCLUSIONS 

In this study we examined the behavior of atoms that 
circulate in a tight toroidal trap. Manipulation of the 
trapping potential along the torus may provide various 
ways to control the stability of the current. The form 
of the potentials that we considered include a step po- 
tential, a sinusoidal potential, a random potential, and 
a gravitational potential. We calculated the value of the 
critical coupling that gives rise to stable persistent cur- 
rents. In addition, we suggested a method that allows the 
detection of circulation. This method involves a probe, 
spatially-dependent potential, which, for different values 
of the circulation gives a different density variation. This 
method requires only a rough measurement of the density 
profile, which may be realizable non-destructively :3(l |. 

Our results indicate that the remarkable progress on 
the manipulation of the trapping potential [37|, and 
more generally the developments in the physics of cold 
atoms may allow us to engineer such mesoscopic systems 
that support persistent currents according to our will. 
Such systems may therefore serve as "superconducting 
switches" with possible important technological applica- 
tions. 
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